Efficiency of quantum and classical transport on graphs 
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We propose a measure to quantify the efficiency of classical and quantum mechanical transport processes on 
graphs. The measure only depends on the density of states (DOS), which contains all the necessary information 
about the graph. For some given (continuous) DOS, the measure shows a power law behavior, where the 
exponent for the quantum transport is twice the exponent of its classical counterpart. For small-world networks, 
however, the measure shows rather a stretched exponential law but still the quantum transport outperforms the 
classical one. Some finite tree-graphs have a few highly degenerate eigenvalues, such that, on the other hand, 
on them the classical transport may be more efficient than the quantum one. 

PACS numbers: 05.60.Gg, 05.60.Cd, 03.67.-a, 



I. INTRODUCTION. 

The transfer of information is the cornerstone of many 
physical, chemical or biological processes. The information 
can be encoded in the mass, charge or energy transported. All 
these transfer processes depend on the underlying structure of 
the system under study. These could be, for example, simple 
crystals, as in solid state physics 1 1], more complex molecu- 
lar aggregates like polymers 0], or general network structures 
[3]. Of course, there exists a panoply of further chemical or 
biological systems which propagate information. 

There are several approaches to model the transport on 
these structures. In (quantum) mechanics, the structure, i.e., 
the potential a particle is moving in, specifies the Hamilto- 
nian of the system, which determines the time evolution. For 
instance, the dynamics of an electron in a simple crystal is de- 
scribed by the Bloch ansatz Q]. Hiickel's molecular-orbital 
theory in quantum chemistry allows to define a Hamiltonian 
for more complex structures, such as molecules |4]. This is 
again related to transport processes in polymers, where the 
connectivity of the polymer plays a fundamental role in its 
dynamical and relaxational properties |5]. There, (classical) 
transport processes can be described by a master equation ap- 
proach with an appropriate (classical) transfer operator which 
determines the temporal evolution of an excitation |2, 6]. 

In all examples listed above, the densities of states (DOS), 
or spectral density, of a given system of size N, 



lie within a semi-circle |8]. As we will show, distinct features 
of the DOS also result in very distinct transport properties. 

II. TRANSPORT ON GRAPHS. 

We start our discussion by considering quantum mechani- 
cal transport processes on discrete structures, in general called 
graphs, which are a collection of N connected nodes. We 
assume that the states \j), associated with a localized exci- 
tation at node j, form an orthonormal basis set and span the 
whole accessible Hilbert space. The time evolution of an ex- 
citation initially placed at node \ j) is determined by the sys- 
tems' Hamiltonian H and reads exp(— «Ht)|j). The classical 
transport can be described by a master equation for the con- 
ditional probability, pkj{t), to find an excitation at time t at 
node k when starting at time at node j. Using also here the 
Dirac notation for a state at node j, the classical time evo- 
lution of this state follows from the transfer matrix T of the 
transport process as exp(Tt) 1 7) . In order to compare the clas- 
sical and the quantum motion, we identify the Hamiltonian 
of the system with the (classical) transfer matrix, H = — T, 
which we will relate later to the (discrete) Laplacian of the 
graph, see e.g. |9, 10]. The classical and quantum mechanical 
transition probabilities to go from the state \j) at time to the 
state |fc) in time t are given by Pk,j{t) = {k\ exp(Tt)|j) and 
TTk,j{t) = |afcj(t)p = I (fc|exp(-iHi)|j)p, respectively. 
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n=l 



contains the essential informations about the system. Here, 
the A„'s are the eigenvalues of the appropriate Hamiltonian H 
or transfer operator T. Depending (mainly) on the topology 
of the system, p{X) shows very distinct features. A classic in 
this respect is the DOS of a random matrix, corresponding to a 
random graph |7]. Wigner has shown that for a (large) matrix 
with (specific) random entries, the eigenvalues of this matrix 
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III. AVERAGED TRANSITION PROBABILITIES. 

Quantum mechanically, a lower bound of the average 
probability to be still or again at the initially excited node, 

7fdiscr(0 = jf ''''ij(^)' obtained for a finite network 

by an eigenstate expansion and using the Cauchy-Schwarz in- 
equality as, 1 1 1], 



TI'dis 



r{t)> 



|adi.cr(i)|'. (1) 



Note that |Sdiscr(i)P depends only on the eigenvalues of H 
but not on the eigenvectors. As we have shown earlier, es- 
pecially the local maxima of Tf{t) are very well reproduced 
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by|a(t)p and for regular networks, the lower bound is exact 
fnll . Therefore, we will use |adiscr(OP the following to 
characterize transport processes. 

Also classically one has a simple expression forp^i^^j.{t) = 

71 Pj.jit)^ see, e.g., U2i, 



Pdu 
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exp I 



Xr,t) 



(2) 



Again, this result depends only on the (discrete) eigenvalue 
spectrum of T but not on the eigenvectors. 

In the continuum limit, Eqs. Q and (|3 can be written as 



W{t) > 



dX p{\) exp(— iA<) 
d\ p{\) exp(— At), 



(3) 
(4) 



The explicit calculation of the integrals is easily done using 
computer algebra systems like MAPLE or MATHEMATICA; 
in many cases, the integrals can also be found in L13J . 



IV. EFFICIENCY MEASURE OF TRANSPORT ON 
GRAPHS. 

Equations allow to define an efficiency measure 

(EM) for the performance of the transport on a graph. We 
stress again, that the EM does not involve any computationally 
expensive calculations of eigenstates. Rather, only the energy 
eigenvalues are needed, which are quite readily obtained by 
diagonalizing H. 

By starting with continuous DOS, since those are mathe- 
matically easier to handle, we define the (classical) EM of 
the graph by the decay of p{t) for large t, where a fast decay 
means that the initial excitation spreads rapidly over the whole 
graph. Quantum mechanically, however, the transition proba- 
bilities fluctuate due to the unitary time evolution. Therefore, 
in most cases also n{t) and |a(t)p fluctuate. Nevertheless, 
the local maxima of reproduce the ones of 7f(t) rather 

well. We use now the temporal scaling of the local maxima of 
|a(f)p as the (quantum) EM and denote the envelope of the 
maxima by env[|a(t)p]. Similar to the classical case, a fast 
decay of env[|a(t)p] corresponds to a rapid spreading of an 
initial excitation. 

For a large variety of graphs the DOS can be written as 



p(A) ^ (AA™ - 



(5) 



with ly > —1 and where Am is the maximal eigenvalue (we 
assumed the minimum eigenvalue to be zero). Since we are 
interested in the large t behavior, p{t) [and also will 
be mainly determined by small A values, such that for t ^ 1 
we can assume p{\) ^ X'^. Than it is easy to show that the 
classical EM scales as 



This scaling argument for long times is well known through- 
out the literature, where 2{1 + i/) = dg is sometimes called 
the spectral or fracton dimension, see, e.g., (3. 

In order to obtain the quantum mechanical scaling for the 
same DOS, we can use the same scaling arguments. For 
t ^ 1, also |a(t)P [and |adiscr(t)P] will be mainly deter- 
mined by the small A values. In fact, for p{X) ^ X" one has 
|a(t)| — p{t). Here, all quantum mechanical oscillations van- 
ish, because we consider only the leading term of the DOS for 
small A. Thus, we furthermore have env[|a(i)p] — \a{t)\'^, 
i.e., the quantum EM reads 



env a 



(7) 



Equation can also be directly obtained from Eq. (|3} with 
p{X) ^ X'^. Of course, Eqs. and Q agree with the so- 
lution for p(t) and env[|S(t)p] obtained from the full DOS 
p(A) ^ (AA„i — X^Y . The same scaling has been obtained 
for the decay of temporal correlations in quantum mechanical 
systems with Cantor spectra 1 15). There, the (full) probability 
7f (t), which was smoothed over time, was used. 

In general, for p{t) ~ 1^''^"^ the exponent Vc\ determines 
the classical EM of the graph because larger Vc\ correspond to 
a faster decay of p{t). Quantum mechanically, we may have 
env[|s(t)p] ~ f~^<i", such that the exponent T'qm determines 
the quantum EM of the graph. Since we consider only the lo- 
cal maxima, the actual (fluctuating) probability 7f (t) [bounded 
from below by might drop well below these values, 

i.e., there ai-e times t at which Tf{t) <^ 1. However, these 
values are very localized in time and the overall performance 
of the quantum transport is best quantified by the scaling of 
env[|aW|2]. 

The difference between the classical and quantum EM is 
given by the factor 



A7'(i)^ln[env[|a(i)n]/ln[p(t)]. 



(8) 



p{t) - 



(6) 



For classical and quantum power law behavior AV{t) is time- 
independent and we have AV = Vqm/Vch Thus, for the 
DOS given above, with v < oo, we get AV — 2, as could 
be expected from the wave-like behavior of the quantum mo- 
tion compared to the normal diffusive behavior of the classical 
motion. 



A. Continuous DOS. 

Two important examples are connected to scaling. An in- 
finite hypercubic lattice in d dimensions has as eigenvalues 
A(ei, . . . , Od) = Et=i A(0«), with A(e„) - 2 - 2 cos e„ 
and Qn G [0, 27r[. Here, one can calculate explicitly 
and Tf{t) and demonstrate that the local maxima really obey 
scaling; we get namely \a{t)\'^ = W{t) ~ \Jn{2t)\'^'^ yiJ. For 
t > 1 this can be approximated by 7f(t) sin^''(2i + 7r/4)/i'' 
fTiil . Since the maximum of the sin-function is 1, the quantum 
measure scales as env[|S(t)p] = cnv[7f](t) ~ t^'^, which is 
what one also obtains from the scaling argument above and 
p(A) ~ A''/^^^. Then — d/2 — 1, and the classical measure 
scales as p{t) ~ t"''/^ (i.e., the spectral dimension is ds = d). 
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As a second example, we take a random graph. It was 
shown that the eigenvalue spectrum of the Laplacian of such a 
graph obeys Wigner's semi-circle law 00], which we obtain 
for 1/ = 1/2 from the DOS given above. For large times both 
measures again obey scaling and we have p(t) ^ t^'^/'^ and 
env[|a(i)p] - t"'^. 

Figure [l] shows the temporal behavior of p{t) and |a(i)p 
as well as the power law behavior of p{t) and env[|a(i)p] for 
(a) an infinite, regular, one-dimensional (ID) graph and (b) a 
random graph. Note that here (and in the following figures, 
too) the very locaHzed minima of do not always show 

up clearly in the logarithmic scale used. 




FIG. 1: (Color online). p{t) and as well as the power laws 

given in Eqs. jTJ and j6j for (a) an infinite regular (ID) graph {v — 
— 1/2) and (b) a random graph whose DOS obeys Wigner's semi- 
circle law {u — 1/2). 

For some DOS, the EMs show no power law behavior 
The DOS given above are bounded from above by a max- 
imal eigenvalue. This does not have to be the case. The 
DOS of small-world networks, for instance, may show long 
A-tails One additional feature of such DOS is that 

they do not obey any simple scaling for small A. Neverthe- 
less, sometimes analytic solutions for, at least, p{t) can be 
obtained LI 6.1 . as, for example for certain ID systems with 
p(A)-A-3/2exp(-l/VA). 

For computational simplicity we consider a 2D system with 

p(A) = A-''cxp(-l/A) (9) 

for A G [0, oo[ and & > 1. The term exp(— 1/A) is usually 
referred to as Lifshits tail, while the term A^** assures that 
limA^co /o(A) = 0. Then, for t ^ 1, the EMs are proportional 
to the product of a stretched exponential and a power law 

|a(t)p =env[|a(t)p] - t(26-3)/2 |^_2V2l) (10) 
p{t) ~ t(2b-3)/4^^p|^_2Vt). (11) 



Furthermore, Tf{t) does not oscillate, an effect which is inter- 
esting in itself but we will not elaborate on this here. 

Although we do not obtain a simple relation between the 
classical and the quantum EMs, [p(t)]'^ and env[|a(t)p] still 
display similar functional forms. Now, however, A'P{t) — 
[2(26-3) lni-8V2i]/[(26-3) is time-dependent. 

Equations (^J and ilQi are only valid for t ^ 1, such 
that limt^oo A7'(i) = V2 for all b. Hence, also here the 
quantum transport outperforms the classical one, which is 
also confirmed by numerical integration of the time-dependent 
Schrodinger equation for a small-world network 1 17]. In fact, 
in both cases the transport is faster than for a regular 2D graph. 
We note that localization is related to other features of the 
DOS as we will recall below. However, at intermediate times 
the quantum EM may drop below the classical EM; the posi- 
tion of the crossover from AV{t) < 1 to A'P{t) > 1 depends 
on the exponent b. 



B. Discrete DOS. 

Up to now, we have only considered continuous DOS, 
where the quantum EM is quicker than the classical one. 
In the following we will consider discrete DOS which are 
obtained by modeling the motion on a given graph classi- 
cally by continuous-time random walks (CTRWs), see, e.g., 
0, and quantum mechanically by continuous-time quantum 
walks (CTQWs) ]9, 10]. The Hamiltonian is given by the (dis- 
crete) Laplacian associated with the graph, i.e., by the func- 
tionality of the nodes and their connectivity. We assume the 
jump rates between all connected pairs of nodes of the graph 
to be equal. 




10° 10^ 10^ 10^ 10^ 
time t 

FIG. 2: (Color online). P(jiscr(*) and |adiscr(i)P for (a) a finite reg- 
ular (ID) graph of size A'^ = 200 with periodic boundary conditions 
and (b) a dendrimer of generation 10 having functionality 2 = 3, i.e. 
N = 3 ■ 2^" — 2. Panel (a) contains also the power-law behavior for 
the infinite regular (ID) graph, see Fig. 0a). 

In general, for finite graphs, Pdiscr(0 ™d env[|adiscr(i)P] 
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do not decay ad infinitum but at some time will remain con- 
stant (classically) or fluctuate about a constant value (quantum 
mechanically). This time is given by the time it takes for the 
CTRW to reach the (equilibrium) equipartitioned probability 
distribution and for the CTQW to fluctuate about a saturation 
value. At intermediate times, Pdiscr(0 ^^'^ env[|adiscr(i)P] 
will show the same scaling as for a system with the corre- 
sponding continuous DOS. Figure |2j a) shows the temporal 
behavior of Pdiscr(0 ™d |adiscr(i)P for a finite regular ID 
graph of size N = 200 with periodic boundary conditions, 
see also j3- At intermediate times, the scaling behavior is 
obviously that of the continuous case shown in Fig.^a). 

Tree-like graphs do not display scaling in general. For 
CTQW on hyperbranched structures (like Cayley trees, den- 
drimers or Husimi cacti), the transition probability between 
two nodes strongly depends on the site j of the initial excita- 
tion . Even in the long time average. 



lim T 



dt TTkjit), 



(12) 



there are transition probabilities which are considerably lower 
than the equipartitioned classical value flollTlll . In Fig.|2jb) 
we display the temporal behavior of Pdiscr(0 ™d |adiscr(i)P 
for a dendrimer of generation 10 having functionality z = 3, 
i.e. N ~ 3 ■ 2^" — 2. Here the classical curve does not show 
scaling at intermediate times. Quantum mechanically, how- 
ever, |adiscr(i)P has a strong dip at short times but then fluc- 
tuates about a finite value which is larger than the classical sat- 
uration value. One should also bear in mind that |adisci (i)P is 
a lower bound and the actual probability will be larger. There- 
fore, according to our measure for intermediate t ^ 1, the 
classical transport outperforms the quantum transport on these 
special, finite graphs. As we proceed to show, the reason for 
this is to be found in the DOS. This is related to (Anderson) 
localization. Anderson showed that for localization the DOS 
has to display a discrete finite series of 5 functions fl^ . 

We consider now a simple star graph, having one core node 
and — 1 nodes directly connected to the core but not to each 
other The eigenvalue spectrum of this star has a very simple 
structure, there are 3 distinct eigenvalues, namely Ai = 0, 
A2 = 1, and A3 — N, having as degeneracies gi — 1, g2 — 
N — 2, and — 1. Therefore we get 



PmscM = ^[l + (^-2)e-* + e-(^-2)* 
7fdi.cr(i) > ^|l + (iV-2)e-"+e-^(^-2)* 



(13) 



.(14) 



Obviously, only the term |(^ - 2) exp(-it)p/iV2 ^ (N - 
2)2/7V^ in Eq. (tlH is of order 0(1). AH the other terms ai-e of 
order 0{l/N) or 0(1/A^^) and, therefore, cause only small 
oscillations (fluctuating terms) about or negligible shifts (con- 
stant terms) from {N - 2)'^/N^. 

Having only one low lying eigenvalue which is highly de- 
generate and no other eigenvalue of a degeneracy of the same 
order of magnitude, results in Pdisci(0 < T^disait) for all 
times t and Pdiscr(^) < |Sdiscr(OP for almost all times t. 



Figure |3] shows the temporal behavior of p^^g^j.{t), 7fdiscr(i), 
and |adiscr(OP for ^ = 10. Now, for all times, the quantum 
transport is slower than the classical one. We also see that 
|adiscr(OP fluctuates about {N - 2 f/N^ = 16/25. 
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FIG. 3: (Color online). Pdiscr(*)' 7i"discr(t), and |adiscr(i)|^ for a star 
with N = 10. 



In general we find for our star-graph that the classical EM is 
lower than the quantum EM. This result is to some extent also 
obeyed by dendrimers and by other hyperbranched structures. 
These, too, have a few highly degenerate eigenvalues, all other 
degeneracies being an order of magnitude less, which results 
in the absence of any scaling of cnv[|adiscr(^)P], see Fig.|2jb). 
Of course, the details are much more complex due to the more 
complex structure, we will elaborate on this elsewhere. 



V. CONCLUSION. 

We have proposed a measure to classify the efficiency of 
classical and quantum mechanical transport processes. De- 
pending on the density of states, the quantum transport out- 
performs the classical transport by means of the speed of the 
spreding of an initial excitation over a given system. For al- 
gebraic DOS, the EMs confirm the difference between clas- 
sical diffusive and quantum mechanical wave-like transport. 
Also for small-world networks the quantum mechanical EM 
is lower than the classical one, i.e. the quantum mechanical 
transport is faster. 

However, for some finite graphs with a few highly degen- 
erate eigenvalues it may happen that the classical transport is 
more efficient, i.e., that the (quantum) states become local- 
ized. We have shown this analytically for a simple star graph. 
More complex structures, like dendrimers or hyperbranched 
fractals, show an analogous behavior. 
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